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Abstract 

An analytical solution for the time evolution of decay of two identical non interacting quantum 
particles seated initially within a potential of finite range is derived using the formalism of resonant 
states. It is shown that the wave function, and hence also the survival and nonescape probabil- 
ities, for factorized symmetric and entangled symmetric/antisymmetric initial states evolve in a 
distinctive form along the exponentially decaying and nonexponential regimes. Our findings show 
the influence of the Pauli exclusion principle on decay. We exemplify our results by solving exactly 
the s-wave 5 shell potential model. 
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I. INTRODUCTION 



Decay is one of the oldest topics in quantum theory. It describes the escape by tunneling 
of particles from a certain region. In addition to the wave function itself, two quantities 
of interest are the survival and the nonescape probabilities. As time evolves, the former is 
defined as the probability that a decaying particle remains in its original state and the latter 
as the probability that the decaying particle remains within the interaction potential. 

In this work, the wave function for the decay of two identical non interacting particles is 
obtained both as a discrete expansion involving the full set of complex poles (and resonant 
states) of the outgoing Green's function to the problem and also as a discrete expansion that 
includes proper complex poles plus an integral contribution which is more appropriate to 
study the long-time be 



single-particle decay 



[i 



ravior. The above approaches have been used into the description of 
2J. There, in general, one finds, after an ultrashort non exponential 
contribution, an exponential decaying regime that depends on a single dominant resonance 
;erm followed at long times by an inverse power law as t~ 3 (s wave symmetrical potentials) 
l|. The short and long-time nonexponential regimes have been experimentally verified in 
recent years [3|, |4|. 

The escape behavior of identical particles is of interest, particularly since recent studies 
have shown that the Pauli exclusion principle has consequences upon the particular power 
law obeyed in the free time evolution (no potential) of identical particles confined initially 
in a finite region of space [5|. The decay of a few-body Tonks-Girardeu gas has also been 
studied in this context by considering an integration over all the coordinates of the particles 
but one. There, it was also found an inverse-power law proportional to t~ 3 [6j. The approach 
followed here considers both the survival and nonescape probabilities as truly multiparticle 
observable quantities and as such the integrations will be all done in one single step. 

The paper is organized as follows. Section [Til reviews succinctly the formalism of resonant 
states, in particular regarding the time evolution of single-particle decay. Section lITTlprovides 
a derivation of exact expressions for the time evolution of decay of two identical particles as 
a resonance expansion in terns of Moshinsky functions, both for factorized symmetric and 
symmetric/antisymmetric states and introduces the two-particle expressions for the survival 
and nonescape probabilities. Section [IV] discusses a model calculation involving the s wave 
S shell potential which involves analytical expressions for the two-particle solutions both 



2 



along the exponential and long-time nonexponential regimes. Finally, Section |V] provides 
the concluding remarks. 



II. FORMALISM 

Consider a single particle confined at t — along the internal region of a real spherically 
symmetrical finite-range potential, i.e., V(r) = for r > a. For simplicity we consider s 
waves and choose as units h = 2m = 1. The solution to the time-dependent Schrodinger 
equation in the radial variable r, as an initial value problem, may be written at time t in 
terms of the retarded Green's function g(r, r'; t) of the problem as 

#(r,t)= ! #(r,r',t)^(r',0)dr', (1) 
Jo 

where \l/(r, 0) stands for the arbitrary state initially confined within the internal interaction 
region. In what follows, for the sake of the simplicity of the discussion, we refer to potentials 
that do not support bound and antibound states. A convenient form of the retarded time- 
dependent Green's function is expressed in terms of the outgoing Green's function G + (r, r'; k) 
of the problem. Both quantities are related by a Laplace transformation where the Bromwich 
contour corresponds to a hyperbolic contour along the first quadrant of the k plane [l| which 
may be evaluated by deforming the integration contour from — oo to oo along the imaginary 
k axis. This allows us to write the retarded Green's function as [l| 

oo 

g(r,r',t)= ^ u p (r)u p (r')M(z p ), (r,r') f <a (2) 

p=— oo 

where the sum extends over the complex poles of the problem, the notation (r,r'Y means 
that the point r = r' = a is excluded in the above expansion (otherwise it diverges) and the 
function M(z r ), the so called Moshinsky function, is defined as [l| 
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-ihk t/2m 

M{z r ) = — I dk = -w{iz r ), (3) 

2vr J_ 00 k - n r 2 

where z r = — exp(— z7r/4)K r .t 1//2 , with r = ±p, and the function u>(C) = exp( — (^ 2 )erfc( — i^) 
stands for the Faddeyeva or complex error function (?}] for which there exist efficient com- 
putational tools [8]. The functions u p (r) in ([2]) correspond to the so called resonant states 
(also known as quasinormal modes) which are solutions to the Schrodinger equation of the 
problem obeying purely outgoing boundary conditions which imply that the corresponding 



energy eigenvalues are complex, i.e., E p = = £ p — where S p yields the resonance 
energy of the decaying fragment and T p stands for the resonance width, which yields the 
lifetime r p = 1/T P of a given resonance level. The lifetime of the system is defined by the 
longest lifetime, i.e., the shortest width. The complex poles ftp — dp Q^-Q 

distributed 

along the third and fourth quadrants of the complex k plane in a well known manner {9]. 

For proper poles, i.e., a p > b p , Eq. fl2]) may be written alternatively, using some properties 
of the Faddeyeva function, as [1] 

00 

g(r, r', t) = V] u p (r)u p (r')e~ i£pt e~^ Vpt + I(r, r'\ t) (4) 
P =i 

where 

00 

J(r, r'; t) = YM p {t)u p {t')M{-z p ) - u;(r) U ;(r')M(z. p }, (5) 
P =i 

and we recall that (r,r'Y < a . Since the potential is real, it follows from time-reversal 
invariance that uJr) = u*(r) and * The last term in (HI) becomes relevant both 

at ultrashort and asymptotic long times [Jj. The description of ultrashort times is more 
involved and will not be considered here since it possibly requires different considerations. 

There is another route to analyze the long-time behavior of g(r,r';t) which follows by 
closing the Bromwich contour mentioned above along a straight line C\ that is 45° off the 
real axis and goes through the origin The resulting expression, that is equivalent to Eq. 
§5§, reads 

I(r, r';t) = (i/n) [ G + (r,r'; k) exp(-ik 2 t)2kdk. (6) 



It turns out that this integral term may be evaluated at long times by the steepest descent 
method as it has a saddle point at k = and hence one may perform a Taylor expansion of 
G + (r, r' , k) around that value to evaluate the integral [l]. Thus, alternatively, at long times, 
the retarded propagator may also be written as 

00 x ^2m-l 

£(2m+l)/2 Q^m- 



* Q2m— 1 

9 (r,r',t) «J3«p(r)«p(0^*e-^+ £ , (2t ^ 1)/2 ^zr ^r, r', k) ^ (7) 



p=l m=l 

where r\\ = 1/ y/Awi, r^ 2 = — a/«/(647t) and r^ 3 = — 1 / V40967U. For decay of a single particle 
it suffices to take m — 1 in Eq. ([7]). However, already for two particles, higher values of m 
are required, as discussed below. 
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III. TWO IDENTICAL PARTICLES 



In the case of a system of identical non interacting particles, it is known that the Hamil- 
tonian H must be symmetric under the permutation of the indices of the particles so the 
exchange operator and H necessarily commute. Thus, it is enough to impose the appropri- 
ate symmetry/antisymmetry on the initial state ^(j/i, 2/2, 0) since symmetry is conserved as 
time evolves. Hence, the time evolution for decay of two identical particles may be written 
as 

#(r,t) = / / g(r 1 ,y 1 ,t)g(r 2 ,y 2l t)^(y,0)dy 1 dy 2 , (8) 



JO JO 

where r and y denote, respectively, (r 1 ,r 2 ) and (yi,y 2 ). 

A simple choice, which corresponds to a symmetric state, is given by the product of single 
particle states ip a {yi) and ip a (y 2 ), with a denoting the state, 

*(y,0)=Mvi)Mv2)- (9) 

Substitution of (jHJ) and <^ into (|Hj) yields the factorized symmetric state 

00 

*(r,t) = C p , a C q , a u p ( ri )u q (r 2 )M(z p )M(z q ) = 

p,q=-oo 

oo \ / 00 \ 

C p , a u p (n)M(z p ) J C q , a u q (r 2 )M(z q ) J (10) 

\p=— oo / \q=— 00 / 

where (r 1; r 2 ) < a and C n>a , with n = p,q, is given by 

/ u n {y)ijj a {y)dy. (11) 







Another choice for the initial state consists of the linear combination of single-particle states 
i[ , s (yi) and ip s (y2)- Here, s = a, (3 refers to the possible states of the two particles, and hence 
we may write 

*(y, 0) = ^=(iJ> a (yi)MV2) ± MvOMv*)), ( 12 ) 

where respectively, the plus sign refers to entangled symmetric and the minus sign to entan- 
gled antisymmetric states. Then, substitution of ( IT2|) and (J2]) into (JHJ) yields 

^M) = -= (Cp^^C^C^u^u^Miz^Miz,), (13) 

" p,q=~oo 
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where (ri,r2) < a and the coefficients C n ^ follow by replacing a for (3 in (II ip . It is worth 
meaning that the coefficients { C7„„ } , which involve only single-particle states, MM, the 
relationship [1[ 

Reffj^J =1, (14) 

where C n>s is defined as ffTTj) with V'sd/) substituted by ip* s {i)). Hence for real initial states, 
C n ,s = C n ,s- Although the C n ^ s are complex and its real part may be negative, they play a 
most relevant role in time-dependent expansions as discussed in Refs. Il0|. 

Alternatively, one may consider Eq. ([7]) instead of ()2]) to calculate Eqs. f lTUj) and f|T3|) . 
This last procedure provides explicit analytical expressions for the exponential decaying and 
long time inverse power terms. We shall consider both possibilities for the model calculation 
below. 



A. Survival and nonescape probabilities 

The survival amplitude of a two-particle system is defined as 



a pa 



A(t)= / **(ri,r 2j 0)*(ri,r2,t)dr 1 dr 2 ; (15) 
Jo Jo 

hence, the survival probability is given by 

S(t) = \A(t)\ 2 . (16) 

The nonescape probability of a two-particle system is defined as 

|*(r,*)| 2 dridr 2 . (17) 

Once ty(r,t) is known, the calculation of S(t) and P(t) follows from Eqs. (I16p and (I17p . 

Notice that P(t) is always larger than S(t). This is a general feature for these quantities, 
that holds also in single-particle decay, which follows from their definition using the Cauchy- 
Schwarz inequality, i.e., in general P(t) > S(t). 



IV. MODEL 

The S shell potential, whose mathematical simplicity allows it to describe the essential 
physical features of the time evolution of decay is suitable for calculations and extends 



the work done for the free-particle case. For two particles this potential may be written as 



V( ri ,r 2 ) = \6{n -a) + \5{r 2 -a). (18) 

For the single-particle case the resonant states of the corresponding problem satisfy the 
Schrodinger equation with complex energy eigenvalues. They read 

{A p sm(K p r) ; r < a 
(19) 
B p exp(iK p r) ; r > a 

and are normalized according to the condition 

/ u 2 p (r)dr + iu 2 p (cl)/2k p = 1. (20) 
Jo 

From the usual boundary conditions for a 5 potential the complex eigenvalues satisfy [l| 

2ik p + \( e 2iK » a - 1) = 0. (21) 

There are well established procedures to calculate the complex poles n n solving fT2~Tj) lj, |2j. 

A convenient feature of this model is that the outgoing Green function G + (r,r', k) may 
be written as the simple analytical expression Q 

exp(ikr ') — (X/k) sin( k(r ' — a) ) exp(z/ca) 



G (r, r ;k) — : — 

k 



(22) 



1 + (X/k) sin(fca) exp(ifca) 

and hence, the partial derivatives of G + (r, r'; k) that appear in Eq. ([7]) may also be written 
explicitly. They are given in the Appendix |A] and are used to obtain the expressions for the 
time-dependent wave solutions discussed below. 

As initial states we choose appropriate combinations of infinite box states 

MV) = \F™(—), (23) 



a v a 

with s = (a, ft). Specifically, for the factorized symmetric state (JSJ) we choose the product 
of two infinite box states with a = 6 whereas for the entangled symmetric/antisymmetric 
states (112"]) we take a = 1 and (3 = 6. 

It is worth mentioning that the state of the system in the Hilbert space corresponds to a 
ket that can be factorized into a space and a spin part as we are assuming a spin-independent 
Hamiltonian. The spin part of the state ket has to be taken into account when considering 
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FIG. 1: (Color online) Plot of the lnS'(i) (lower solid line) and lnP(t) (upper solid line) as a 
function of time in lifetime units for the factorized symmetric state using Eq. (|10[) . The inset shows 
a calculation of the same quantities at short times using the purely exponential contributions (solid 
lines) and its corresponding comparison with exact calculations in terms of Moshinsky functions 
(dotted lines). See text. 

the total symmetry of the state. In this work we will address only the spatial part of the 
state ket. We shall consider calculations involving both the exact analytical and asymptotic 
expressions for the retarded propagator. 

The calculations below are made using both the exact expressions ( fTUl) and (|T3l) and the 
corresponding approximate expressions that exhibit the long time behavior explicitly. In the 
next subsection we provide explicit analytical expressions for the corresponding symmetric 
and antisymmetric wave functions. 

The potential parameters employed in the calculations are: A = 6 and a = 1 and it was 
sufficient to consider 20 poles. 

A. Symmetric space wave functions 

We first consider the factorized symmetric wave function given by Eq. ( fTOj) . It follows 
that the long time contribution of g(r, r', t), given by the second term on the right-hand side 
of (J71), may be truncated at order t~ 3 ^ 2 and thus the wave function reads 



S 



*(r,f) « ^C p ^(7 fta ti p (r 1 )« a (r 2 ) e ^ + «')* e -it r * +r «) t 



(1 + Aa) 4 t3 



" /L ^^^(^(n) +r 1 u p (r 2 ))e-^V^, (24) 



(1 + Aa) 2 t 3 / 2 

where (ri,r 2 ) < a and we have used the shorthand D a = f yip a (y) dy. Equation ()24l) is a 
description up to leading terms of the solution; it has essentially three contributions: purely 
exponential decaying terms, a purely non-exponential inverse power law and a mixed term, 
made up by several inverse-power multiplied by exponential contributions. Figured] exhibits 
a plot of both lnS'(t) and lnP(t) in units of the lifetime of the system, (which is given by 
7~i = 1/Ti, since 1^ is the shortest decay width). We observe several regions of interest for 
both the survival and nonescape probabilities: for short times there is an exponential regime 
in which the slope of lnS'(t) and lnP(t) goes as — 2r 6 , then, there is a change of the slopes 
into — 2T\ (since this is the exponential term with the longest duration); this is followed 
by interference contributions between the purely exponential terms and the purely inverse 
power term and finally the non-exponential asymptotic regime that goes as t~ 6 . 

We now consider the entangled symmetric wave function (plus sign) given by Eq. ffTSj) . 
It follows that the long time contribution of g(r,r',t), given by the second term on the 
right-hand side of ([7]), may be truncated at order t~ 3//2 and thus the wave function reads 



V p,q=l 



( D pC P , a + D a C p ,p)(r 2 u p (rt) + r 1% (r 2 )) e -^e-^, (25) 



v/2(i + Ao)W jr{ 

where we recall that (r 1; r 2 ) < a. We note that Eq. (125!) is again a description up to leading 
terms of the solution; it has essentially three contributions: purely exponential decaying 
terms, a purely non-exponential inverse power law and a mixed term, made up by several 
inverse-power multiplied by exponential contributions. In Fig. [2] we see several regions 
of interest for both lnS'(t) and lnP(t): for short times there is an exponential regime in 
which the slope of lnS'(t) and lnP(t) goes as —(Fx + T 6 ) rs — r 6 (since T e ^> Fi), then, 
there is an interference among the purely exponential terms and between them and the 



mixed term; the next region goes with slope — 2Ti; then there is an interference between the 
purely exponential terms and the purely inverse power term and finally the non-exponential 
asymptotic regime that goes as t~ 6 . 



o.o 




FIG. 2: (Color online) Plot of the In S(t) (lower solid line) and lnP(t) (upper solid line) as a 
function of time in lifetime units for the entangled symmetric state using Eq. (|25p . The inset shows 
a calculation at short times of the same quantities using the purely exponential contributions (solid 
lines) and the corresponding comparison with exact calculations in terms of Moshinsky functions 
(dotted lines). See text. 

B. Antisymmetric space wave function 

The entangled antisymmetric wave function follows by choosing the minus sign in (I12p 
and proceeding in a similar fashion as above. At long times, the propagator g(r,r',t) has 
to be expanded up to t~ 7 / 2 order since the terms of order t~ 3 and t -4 in \I/(r, t) cancel out 
exactly. Thus, we obtain a time-development given by 
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1 oo 
P,9=l 



q,a)Up(.ri)U q (r 2 



]e -i{s p +e q )t e _i(r p +r,)t + 



2 IO771773 \ {rfr 2 - r|n) (DpG a - G P D C 



y/2(l + Xa)H 5 

00 

v^(l + Aa)W» E^^" - A*C^)(r 2 u p (n) - nu^e^e'^ (26) 
where (ri,r 2 ) < a and G Q = J^y 3 ip a (y) dy. Note that Eq. ( |26|) has also three leading 




FIG. 3: (Color online) Plot of the In S(t) (lower solid line) and In P(t) (upper solid line) as a function 
of time in lifetime units for the entangled antisymmetric state using Eq. (|26jl . The inset shows a 
calculation of the same quantities at short times using the purely exponential contributions (solid 
lines) and the corresponding comparison with exact calculations in terms of Moshinsky functions 
(dotted lines). See text. 

contributions: a purely exponential decaying one, a purely non-exponential inverse power law 
(now of order t~ 5 ) and again a mixed term. It is relevant to realize that {C p>a C q ^ — C p ^C qtQ ) 
enforces an exact cancelation of the terms in the sum whenever p = q and hence there are 
not contributions of order — 2T p , for any p, in the exponential decaying terms as occurs in 
the symmetric case. In Fig. [3] we can identify several characteristic regions for lnS'(i) and 
In P(t): for short times there is an exponential regime in which the slope of both hiS(t) 
and lnP(t) goes as — (ri + T 6 ) m — T 6 (recalling that T 6 3> Ti); next one observes a short 
linear region of slope — (fi + T 2 ); then it follows an interference term between the purely 
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exponential and the mixed contributions; the next region comes from the mixed terms and 
goes with slope — Ti\ then there is an interference contribution between the mixed terms 
and the inverse power term, and finally the non-exponential regime as t~ 10 . 

A salient feature of the insets in Figs. [U [2] and [3] is that the survival probabilities are 
much less than the corresponding nonescape probabilities. This also occurs for single-particle 
decay {2]. It occurs whenever the decay process starts from an excited state. Otherwise along 
the exponentially decaying regime S(t) and P(t) are practically indistinguishable. 

In the absence of a potential the problem reduces to the free time evolution of an initially 
confined two particle state. This has been recently discussed in Ref. jlj]. Their results are 
reproduced by letting A = in Eqs. (l25l) and ([25]), in which case there are no pole contri- 
butions and the symmetrical and antisymmetrical free wave solutions evolve respectively as 
i~ 3 and t~ 5 . 



V. CONCLUDING REMARKS 



This work shows that the character of the time evolution of decay of two identical non- 
interacting particles depends on whether the initial state is factorized symmetric or entangled 
symmetric/antisymmetric. The results obtained indicate that each of the above initial cases 
exhibits in general a distinctive behavior along the exponential and long-time nonexponential 
regimes. The differences in the exponential decaying regime may be observed in the different 
slopes of the plots discussed here and, for the asymptotic long-time regime, in the different 
inverse power laws obeyed. 

We have restricted the discussion to expansions of the time- dependent solution along 
the internal interaction region. This is sufficient to calculate the survival and nonescape 
probabilities. One might also consider the time evolution for decay along the external 
interaction region by generalizing to two particles the resonant expansion discussed in Ref. 

which yields for the single-particle time evolving retarded Green's function the expression 

00 

g(r,r',t)= y~] u p (a)u p (r')M(z p ), r' < a, r > a, 
p=— 00 

where the argument z p of the Moshinsky function now reads z p = exp(-?7r/4)(l/2t 1 / 2 [((r - 
a) — 2n p t]. On the other hand, one may extend the discussion given here to include the 
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effect of bound states just by adding the corresponding bound and antibound terms to the 



resonance sums ffTUj) and (fTB"]) jl]. This may lead to interesting effects as 'trapping' |12|. 

As our example exhibits, the survival and nonescape probabilities both go at asymptot- 
ically long times, respectively, as t~ 6 for both factorized/ 'entangled symmetric states, and 
as t~ 10 for entangled antisymmetric states. In the limit of a vanishing interaction potential 
our results tend to the free evolving case discussed in Ref. [5(. The above long-time re- 
sults are in contrast with the t~ 3 behavior for decay of a single particle [1]. Evidently, the 
exponential decaying regime is the one more accessible to experimental verification. Here 
the factorized symmetric states possess a different behavior than for the entangled symmet- 
ric/antisymmetric states, which, on the other hand, exhibit at short times a similar decay 
rate. 

We hope that the distinct behaviors for the time evolution of decay discussed here might 
be experimentally verified in quantum systems where one may manipulate to certain degree 



the potential parameters such as in cold atoms 



a. 



Note- After submission of this work, we became aware of Ref. 13] which develops an 



approach to many particle decay addressing only the long-time decay regime. 
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Appendix A: Derivatives of G + (r,r';k) with respect to k evaluated at k = 

We write down below the explicit expressions for the derivatives with respect to k of 
G + (r, r'\ k) at k = which are required to obtain the analytical expressions for the time- 
dependent solutions given in Subsecs. IIV Al and IIVBI Here we recall that (r, r'y < a. 



d „ , , , , s irr' 



dk G+{ry ' k) 



G + (r, r',k) 



k=0 hi 

irr'(hi(r 2 + (r') 2 ) — h 2 ) 
k=o h\ 
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Q5 



G + (r,r',k) 



where 



k=0 



irr 
3hJ 



(/i 2 (3r 4 + 3(r') 4 + 10r 2 (r') 2 ) + /i 3 (r 2 + (r') 2 ) + fc 4 ) 



/ii = (1 + Aa) 5 



/i 2 = 8Aa 3 + 2A 2 a 4 



-120A 3 a 5 - 80Aa 3 - 20A 4 a 6 - 180AV 



h A = 192AV - 96Aa b + 24AV + 432AV 
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